Astronomy & Astrophysics manuscript no. muratorio 


5th February 2008 


(DOI: will be inserted by hand later) 





Analysis of the circumstellar environment of the B[e] star 

HD 45677 (FS CMa)* 

G. Muratorio 1 , C. Rossi 2 , M. Friedjung 3 

1 Laboratoire d'Astrophysique de Marseille, 2 Place le Verrier, FR-13248 Marseille Cedex 04, France 
e-mail: gerard.muratorio@oamp.fr 

2 Dipartimento di Fisica , Universita "La Sapienza", Piazzzale A. Moro 2 00185, Roma, Italy 
e-mail: corinne .rossi@romal . infn. it 

3 Institut d'Astrophysique, 98 bis Boulevard Arago, FR-75014 Paris, France 
e-mail: fried@iap.fr 

the date of receipt and acceptance should be inserted later 

Abstract We studied the circumstellar environment of the B[e] star HD 45677 through the analysis of the emission lines from 
ionized metals. We used the statistical approach of the self absorption curve method (SAC) to derive physical parameters of the 
line emitting region. Fhe Fe II and Cr II double-peaked emission line structure is explained by the presence of a thin absorption 
component red shifted by ~ 3 kms~'. Fhis absorption component can be interpreted geometricaly as being due to infalling 
material perpendicularly to the disk seen nearly pole-on, as indicated by the emission line structure. Fhe Cr II and Fe II emission 
lines have a complex structure with two (narrow and broad) components, of 45 and 180 km s -1 FWHM for the permitted lines 
and 25 and 100 km s" 1 FWHM for the forbidden ones, respectively. From our best data set of 1999, we obtained a Boltzmann- 
type population law whose excitation temperature is 3 900^'jJJJ K and 3 150^^ K for the narrow component lower and upper 
levels respectively. We obtained an excitation temperature of 3 400*HS K for the broad component upper levels. Fhe forbidden 
lines are found to be formed in the outer regions with higher excitation temperatures of 10 500 ± 1 000 A* and 8 000 ± 1 500 A" for 
the narrow and broad components respectively in 1999. Our results are consistent with line formation in a rotating disk, around 
a young star. In the framework of a very simplified geometrical model, we argue that the narrow components are principaly 
emitted by an optically thin disk seen nearly pole-on, in a region whose minimum radius is estimated to be 4 10 12 c/n, while the 
broad ones are formed in a disk-linked wind. 
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1. Introduction 

HD 45677 (FS CMa) is generally recognized as an early type 
star, usually classified as B[e], or more precisely as showing the 
B[e] phenomenon, even if it is difficult to determine its exact 
evolutionary stage (de Winter et al. 1997] Lamers et al ll998> . 

Fhe presence of a disk, and even more of a bipolar envi- 
ronment, was inferred by polarization measurements (Coyne 
& Vrba lT976l Schulte-Ladbeck et al. lT992l and by UV spec- 
troscopic measurements of accreting gas (Grady et al. I1993> . 

Send offprint requests to: G. Muratorio e-mail: 
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* Based on observations obtained at the Haute Provence 
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These last authors agreed on the presence, around the star, 
of an actively accreting circumstellar, protoplanetary disk, and 
the presence of low-velocity absorption profiles in Fe n lines, 
among other species, while optical spectroscopic measure- 
ments of strong absorption cores of a few lines are commented 
on by de Winter et al. J1997i . 

In the optical HD 45677 is characterised by a rich emis- 
sion line spectrum, including the Balmer lines of hydrogen, and 
both permitted and forbidden lines of neutral and ionized met- 
als (Mg i, Mg ii, Na i, Cr n, Si n, Mn n, Ti n, [Ni n], [S ii] and 
especially Fe ii and [Fe n]). 

Israelian et al. (1996) found a surface gravity of log g = 
3.85 + 0.05 (luminosity class V) from the photospheric wings 
of Hy and HS, while examination of the He i absorption lines 
suggested a spectral class of B2 (T eff =22 000 + 1 500 K). The 
effective temperature was less certain than the surface gravity, 
due to a circumstellar contribution to the He i lines. Israelian 
et al. J1996b also found that the photospheric Si n 4128 and 
4130 A lines were not blended, and could be only fitted by 
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models with a projected rotational velocity v sin ; of less than 
TOkms -1 . 

Cidale et al. J200H used the method of Barbier, 
Chalonge & Divan, based on the magnitude and position 
of the Balmer photospheric (not circumstellar) discontinuity. 
They obtained a spectral type of B2IV-V, corresponding to 
T eff =21 500 ± 300 K and log g = 3.89, in agreement with the 
previous determinations. 

De Winter & van den Ancker J1997i studied the behaviour 
of HD 45677 during the last 25 years, and concluded that the 
large photometric variations of the star are best explained by a 
variable obscuration from circumstellar grains. They also de- 
rived a stellar instrinsic visual-ultraviolet energy distribution 
close to that of a B2 V star, with an infrared excess due to the 
dust emission. The recent spectral variations were discussed by 
Israelian et al. (1996i and Israelian & Musaev ( 119971 . 

In the present study, we analyse the emission line spectrum 
in order to get an insight into the physical processes of line ex- 
citation and geometrical conditions of their forming regions. It 
is clearly impossible at the present stage to calculate a detailed 
synthetic spectrum for this badly known object, so in order to 
analyse the spectrum, we use the statistical approach offered 
by the Self- Absorption Curve method (Friedjung & Muratorio 
[T9"87l Muratorio & Friediung [T988l Baratta et al. [T998l Kotnik- 
Karuzaet al. l2002l 

2. Observations and Reductions 

The spectra of HD 45677 were obtained in various contiguous 
domains in the course of different observing runs at the Haute 
Provence Observatory (OHP) by means of the Aurelie spectro- 
graph attached to the Coude focus of the 1.52 m telescope. 

Until 1999 Aurelie used as receiver a Thomson TH 
7832 photocell "Double Barette", and since 2000 a new 
CCD/EEV which improved the instrument performances (see 
http://. www.obs-hp.fr/www/guide/aurelie/annonce.htmll. 

We obtained many spectra in selected regions; A higher res- 
olution was used in 1999 and 2001-2002 in order to charac- 
terise better the profile of the metallic spectral lines. A detailed 
log of the observations is given in Table. 1 1 (online material). 

The 1995 observations were used in 1996 for the analysis 
of the He i lines 4387, 4471, 5876 and 4438 A (Israelian et al. 
19961. 

The spectra were reduced at the Marseille Observatory us- 
ing procedures runing in the MIDAS environment. 

After the files were calibrated in wavelength, and reduced 
to the heliocentric standard of rest, we identified the spec- 
tral features present in all our spectra. No large variations are 
present in the metallic emission lines within our observing time 
target (see Fig.0. Higher resolution spectra of HD 45677 were 
obtained from October 14, 19 and 20 in 1995 by Musaev with 
the Coude Echelle Spectrometer at the 1-m telescope of the 
Special Astrophysical Observatory of the Russian Academy of 
Sciences (Israelian & Musaev Tl997> . and kindly provided for 
us by Garik Israelian. 

For the purpose of our analysis of the emission lines we 
need the absolute values of the fluxes. To do this, we have nor- 
malized the individual spectra to the local continuum level and, 



following de Winter et al. J1997I . we multiplied the line fluxes 
by the energy distribution of the Kurucz ( fT979l[T993l model of 
T e ff = 22 500K and log g = 4.0 with solar abundance, normal- 
ized to the visual magnitude of the star, Vq=6.6. The results 
are not dependent on the model choosen in the limit of the pre- 
cision of the determination of the star's type by various authors 
(see Introduction). We adopted the same continuum distribu- 
tion for the three studied epochs. 

2.1. Line identification 

Extensive line identification lists of HD 45677 are available 
from the literature (Swings I1973I Andrillat et al. 1 19971 de 
Winter et al. 119971 1. However, they generally refer either to 
lower resolution spectra and/or to limited spectral ranges. We 
therefore prepared a new line list, based on our observations, 
which is presented in Table. 12 (online material). Emission and 
absorption components identified in different species are sum- 
marized in Table. 13 (online material). 

2.2. Line components 

Figure shows that the lines of the most numerous species, 
namely Cr n , Fe n and [Fe n ], did not significantly vary dur- 
ing 1995-2002. This will be best illustrated at the end of this 
section. 

As shown in Figure^ and, more in detail, in Figures|3|3El 
and 17 (online material), in the higher resolution spectra of HD 
45677, the strongest permitted lines present a double peaked 
emission. This feature has been in previous papers ascribed to 
emission from a rotating equatorial disk (de Winter et al. H997i . 
But they can be as well attributed to absorption in the emitting 
medium or in a circumstellar shell according to Grady et al. 
dl 993i > from their ultraviolet observations. 

Israelian and Musaev (I1997> discussed the narrow absorp- 
tion components observed in Mg n A 448 1 A and in the hy- 
drogen lines in the framework of both infalling material with a 
patchy structure, and outfall-wind. They also noticed the strong 
variations of the Mg n line from pure absorption to P Cygni 
profile on timescales of a few days. 

Actually, in the higher resolution spectra the central dips of 
the strong Mg n A 4481 A and Cr n A 4558.6 A and A 4588.2 A 
go well below the continuum level (Figs.[0and|2j. 

In our lower resolution 1995 spectra (top), the absorption 
is probably smoothed out, but is nevertheless present, as con- 
firmed by the fitting procedure discussed below. The other 
emission lines are either sligthly flat-topped or present a dip at 
the red side of their maximum, again suggesting the presence 
of an unresolved central (or slightly red-shifted) absorption. 

Indeed, the line shape is that typical of an emission line 
with a central absorption, rather than that of a rotating disk seen 
close to edge-on and we shall see in Section|4]that these profiles 
can be interpreted as originating in a nearly face-on seen disk. 

At the same time, the strongest emission lines clearly 
present extended wings, that, however, are difficult to detect 
in many cases because of line blending. 
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Figure 1: The spectrum of HD 45677 during 1995-2001. The spectrograms, normalized to the continuum and vertically shifted, 
are, from top to bottom, of 14/01/1995 (R=7000), 19/10/1995 (R=36000), 2/2/1999 (R=15000) and 5/12/2001 (R= 15000). All 
the spectra were obtained at OHP with AURELIE, except the October 1995 one which is from the Russian Special Astrophysical 
Observatory. On the higher resolution spectrograms the deep Cr n absorptions at 4558 and 4588 A show up clearly. 



In Fig.|5]we show in more detail the two most intense lines 
of Cr ii multiplet 44. The line profiles are well reproduced with 
a multiple Gaussian fit, with two components in emission with 
a FWHM of 45 and 180 km s _1 FWHM respectively, having the 
same central wavelength, and one absorption with a FWHM of 
25-30 kms" 1 , slightly displaced to the red. (See the infalling 
material model discussed by Israelian and Musaev.) 

Strong absorption components have also been observed in 
the Na i D doublet since a long time, and were discussed by de 
Winter et al. J1997i . In our spectra the Na i D doublet shown 
in Fig. 17 underwent dramatic variation in both intensity and 
position of the emission components, as well as in the intensity 
of the central absorption (its position is stable within the fitting 
precision). This behaviour is not unexpected because Israelian 



et al. |T996| report variability for these lines within 2-3 days. 
Due to the intense absorption, the narrow emission component, 
if present, cannot be found by the fitting procedure. The ve- 
locity and fluxes of absorption and broad emission components 
are given for the three epochs in Table. [2 



Table 1: Flux (in 10" 13 ergcm 2 s ) and velocity (in kms ') 
variations of the Na I Dl and D2 components. 





1995 
flux 


1995 
vel 


1999 
flux 


1999 
vel 


2002 2002 
flux vel 


NalDl a 


-0.88 


31 


-1.85 


24 


-2.60 24 


NalDl e 


2.77 


48 


4.84 


4 


2.86 9 


Na I D2 a 


-1.26 


31 


-2.05 


22 


-2.55 22 


Na I D2 e 


3.34 


48 


5.24 


5 


2.88 6 



As for the Fe n lines, a central absorption is clearly present 
in the line profiles at A 4508.288 A in Fig. □ and Fig. and 



at A 4520.224 A, 4549.474 A and 4583.837 A in Fig.ffl The 
absorptions, like the Cr n ones, are red shifted an average of 3 
km s~' with respect to the central wavelength of the emission; 
the shift is the same for all the observing runs (Table. 

The same set of components than for Cr n lines with similar 
parameters is present in the Fe n lines (FWHM of 45 and 180 
kms -1 ), while the [Fe n] display narrow and broad emissions 
with smaller FWHM of 25 and 100 kms 1 , respectively (see 
Fig. for highly blended lines and Fig. 0] for more separated 
ones). 

As the FWHM of the Fe n absorption components FWHM 
are similar to the Cr n ones (^30 kms -1 ), these components 
probably originate from the same region. 

We postulate that in HD45677 the double-peaked Fe n line 
structure reported in the literature is in fact due to the same type 
of absorption and originates from material in the star's line of 
sight. Their redward displacement can be interpreted as due to 
the infall of the absorbing material. This result is in agreement 
with the observations and the conclusions of Grady et al. (I1993> 
and Israelian and Musaev (119971 . 

The two emission components clearly originate from re- 
gions of different dynamical behaviour, a possible model of the 
environment of FS CMa will be given in Sec.0] 

Based on the above results, we have fitted all, strong and 
weak, metallic permitted lines identified in our spectra with 
three Gaussian line components, central absorption narrow 
emission and broad emission. The forbidden lines were fitted 
with two emission components. The same central wavelength 
was assumed for the two emission components. The derived 
line fluxes, FWHMs and central radial velocities are given in 
Table. 12 (online material). For the line fluxes we estimated 
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Figure 2: The two most intense lines of multiplet 44 of Cr n 
observed the 5/12/2001, showing the central absorption well 
below the continuum and the three Gaussian components' fit 
. Each line profile (solid line) is fitted (dotted line) by three 
Gaussian components (dashed line). The residual (long dashed 
line) appears around the continuum (straight line). Notice the 
rising of the residual towards shorter wavelengths which is in 
both cases due to the presence of a strong Fe n line. 

from the fitting procedure an error from 15% for the strongest 
lines to 50% for the weakest lines. 



Table 2: Mean radial velocities (in km s _1 ) of the most numer- 
ous species, absorption and emission of Cr n, Fe n, and [Fe n] 
emission. The number of lines is given in parentheses. 





1995 




1999 




2001 




Cr II a 


25.4 ± 


1.3(10) 


25.4 ± 


1.4(11) 


21.8: 


t0.8(4) 


Cr II e 


23.6 ± 


1.3(11) 


23.1 ± 


1.2(11) 


18.6 


+ 2.0(4) 


Fe II a 


29.4 ± 


0.7(58) 


28.2 ± 


0.6(59) 


24.2 


+ 0.6(35) 


Fe II e 


25.8 ± 


0.5(69) 


24.5 ± 


0.4(64) 


21.6 


+ 0.5(47) 


[Fe II] 


21.5 ± 


0.5(85) 


22.5 + 


0.4(77) 


21.4 


+ 0.3(32) 
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Figure 3: Profiles and various components of the blend 
4508.288 A (Fe n multiplet 38), 4509.610 A ([Fe u] multiplet 
6F) on the 2/2/1999 (see Fig.|2]for the explanation of the differ- 
ent curves). The different component FWHM are given in the 
text. 
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Figure 4: Profiles and various components of 5273.38 A ([Fe n] 
multiplet 18F), 5276.002 A (Fe n multiplet 49), and 5278.374 A 
([Fe n] multiplet 35F) on the 3/2/1999 (see Fig.|2for the expla- 
nation of the different curves). The different component FWHM 
are given in the text. 



The mean heliocentric radial velocity of the narrow emis- 
sion and central absorption components of Cm, Fe n and [Fe n] 
in the three epochs are given in Table|2] 

We compared the fluxes of the narrow emission compo- 
nents of Fe ii and [Fe n] lines in 1995 and in 2001-2002 with 
those of 1999. 

Between 1995 and 2001-2002 the fluxes of the permitted 
lines are the same within the measurements errors, except for 
the weaker lines which are fainter in 1995. In December 2001- 



January 2002 the forbidden lines have weakened by somewhat 
less than 30% with respect to 1995 and 1999. As for the broad 
components, their larger measurement errors only allow to put 
an upper limit of 30% to their variation during our observing 
period. 

In the following we shall analyse the broad and narrow 
emission components of these species, with special regard to 
the February 1999 data which cover a wider spectral range. 
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3. Analysis of the emission line fluxes 

The spectrum of HD 45677 is rich in emission lines (see on- 
line material table 13), whose strength can be used to derive 
information about the physical parameters and geometry of the 
circumstellar region. Here we analyse the intensity of the per- 
mitted and forbidden ionized iron emission lines and permitted 
ionized chromium emission lines in a statistical way, the Self- 
Absorption-Curve Method (SAC), following the procedure de- 
scribed by Friedjung & Muratorio d!987> . and used among oth- 
ers by Muratorio et al. dl992> . Viotti & al. J1999L Viotti et 
al. (|2000|, Kotnik-Karuza et al. (2002 1 and van den Ancker et 
al. ( EUOH . Details on the SAC method are given in the SAC 
Manual (Baratta et al. 119981 available at www.rm.iasf.cnr.it 
/ftp/uvspace/). 

The Self absorption Curve Method has been discussed by 
Kastner ( 1999 1, in the framework of the escape probability for- 
malism of radiative transfer. He found, with such a formalism, 
that for a static source having negligible differential velocities, 
the method is valid up to optical thicknesses of about 10. The 
self absorption curve we use is very similar to that found from 
Kastner's approach. As we shall see, the use of a static model 
can be justified for one of the emission line components, while 
most of the Fe n emission lines examined, do not have an opti- 
cal thickness larger than 10. 

The method makes a certain number of assumptions, as fol- 
lows : 

(a) The populations of different levels belonging to the 
same spectroscopic term are approximately proportional to 
their statistical weights. 

(b) The self absorption of lines of different multiplets is 
always the same function of the optical thickness. In the opti- 
cal region studied by us, our analysis is mainly based on lines 
whose lower levels are both even and metastable, while the up- 
per levels are odd. The processes of excitation are supposed to 
be similar. If such a situation exists, we may expect that effects 
of inhomogeneities in the emitting region, affect the different 
lines in almost the same way. 

(c) Selective excitation mechanisms, and particularly 
pumping of populations of upper levels of emision lines be- 
cause of wavelength coincidences with other strong emission 
lines, are fairly infrequent and can be neglected for most lines. 
Such an assumption would be quite dangerous in the ultraviolet 
wavength range (Eriksson et al. l2004> . 

Other assumptions can often be made, as in this paper, 
when applying the method: 

(d) The populations of different collisionally excited even 
lower terms of Fe + follow roughly Boltzmann's law. In fact one 
can expect LTE at least for those having a lower excitation po- 



tential. The only relevant calculations of Verner et al. Q002 ) 
for emission by the B+D Weigelt blobs of r\ Car, with an es- 
timated electron density of the order of only 10 6 cm~ 3 , suggest 
that Boltzmann's law is obeyed up to 2.3 ev, with deviations as 
a function of excitation potential for higher metastable levels. 
Our lines are formed in a circumstellar region, with a presum- 
ably much higher electron density and larger collisional excita- 
tions and de-excitation rates of the metastable levels, as indeed 
is consistent with the determined column densities and limits 



on the size of the line emitting region, given below (see Section 
5.2, for what may happen, if this assumption is relaxed.). 

(e) The populations of the upper odd levels are not in LTE, 
but we suppose that they tend to cluster around a Boltzmann 
law like distribution. Our analysis of HD 45677 is consistent 
with such an assumption for levels, which are not highly ex- 
cited (see below), irrespective of what the dominant excitation 



mechanism may be. Let us note that the Verner et al (|2002i 
i] Car analysis, indicates departure coefficients, which tend to 
be mostly a function of the excitation potential. 

3.1. Fe ii permitted emission lines 
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Figure 5: Upper panel : Log normalized fluxes versus log gfA 
plots for the narrow components of Fe n multiplets 27 (circle), 
38 (square),48 (triangle) with common upper term z 4 D° in the 
1999 spectrum of HD 45677. Lower panel : log normalized 
fluxes versus log gfA - p^i plots of multiplets 27, 38 and 48. 
The multiplets are horizontally shifted according to the popu- 
lation ratio of their lower terms by p = 1.3. 

We first consider the narrow emission components of Fe n. 
For the analysis we adopt Kurucz' (119931 oscillator strengths. 
In order to look for self-absorption effects, we have analysed 
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the narrow component line fluxes of individual Fe n multiplets 
in the [log gf/i, log 24- ] diagram, where the quantity E4- is the 
normalized flux of an individual emission line, and gfA is pro- 
portional, for a given multiplet, to the optical thickness at the 
line centre (see Viotti et al. $20001). This procedure requires of 
course the first of the above assumptions (a). 

In the upper panel of Fig.[3]we show the diagram for mul- 
tiplets 27, 38 and 48 that have z 4 D° as common upper term. 
In the diagram, if the lines were optically thin, the normalized 
emission line flux should be constant within a multiplet. We 
note on the contrary that the strongest lines of multiplets 27 
and 38 have a normalized flux lower than that of the weaker 
multiplet lines. This indicates that they are affected by a certain 
amount of self-absorption. It should also be noticed that we 
have not found for any multiplet any systematic trend, within 
the measurement errors of the line fluxes, with respect to the 
quantum number J, which would suggest a contradiction with 
assumption (a), about level populations within the upper and 
lower terms of each multiplet. Anyhow, a thourough study of 
this effect will be difficult for such complex line profiles, even 
with much higher quality material. 

The bend of these diagrams helps to determine the relative 
population of the lower terms of the multiplets, by their hori- 
zontal shifts with respect to each other by an amount equal to 
the logarithm of this population ratio (e.g. Kotnik-Karuza et al. 
120021 . This has been performed in the [log gfA - pxu log ^-] 
diagram (lower part of Fig. 0» where the multiplets have been 
overlapped assuming a power law with exponent -p. A value of 
p =1.3 + 0.4 eV~ l ) is derived for the exponent of the power law 
of the lower terms population of the multiplets 27, 38 and 48. 

Similarly, by comparing the 1999 diagrams for multiplets 
28, 37 and 49, having z 4 F° as common upper term, we derive 
a mean lower level population gradient of p = 1 .0 + 0.4 eV~' . 

In general these results are fairly uncertain because of the 
small energy range of the lower terms, and the same is the case 
of multiplets with common lower term, which cannot be prop- 
erly used to derive the mean excitation of the upper terms be- 
cause of the small energy range of our data. In order to over- 
come this problem, we have assumed that all the multiplets fit a 
common theoretical self-absorption curve Q (t), and have used 
the least squares method to determine the parameters p, X c , q, 
Y c , needed to overlap the lines of all multiplets to the theoret- 
ical curve in an [X - X c , Y - Y c ] diagram, where they corre- 
spond to the following physical definitions : 

X = loggfA-pCtf-ATc) (1) 

Y = log — + q(Xu ~Xc) (2) 
gf 

Xc = 1.576- log — + log v c (3) 

gc 

v . S ' i Nc . 27re2h 

Y c = log — + log — + w - log (4) 

d 2 g c m e 

and log - 16.977 in cgs units. 



The quantities^/ and^„ are the excitation potentials (in eV) 
of the lower and upper levels of each transition, % c is a char- 
acteristic excitation potential taken different from eV in the 
case of Fe n permitted lines (and equal to 2.75 eV) in order to 
avoid problems of the populations of the higher metastable lev- 
els, whose possible large departures from LTE are pointed out 
by Verner et al. (I2002> . The unknown parameters p and q de- 
scribe the mean population distributions of the lower and upper 
terms, respectively (Muratorio et al. I1992> . In this procedure 
we assume the first of the above assumtions (a), as verified in 
each multiplet [log gfA, log £4-] diagram, and that all the lower 
and upper levels follow a Boltzmann-type law with excitation 
temperatures Ti and T u of 5040/p and 5040/q, respectively. 

The zero points X c and Y c , are related to the Fe + column 
density -j- at the arbitrary choosen level Xc, taken to be equal 
to 2.75 eV for Fe n permitted lines, and to the projected ex- 
tension of the emitting volume S' (S' = S cos i), through the 
equations and (Friedjung & Muratorio 1987 Baratta et al. 
119981 . The column density — also depends on the parameter 
v c (eq. @, which is the velocity broadening of the opacity pro- 
file (see Section l5~2l . For the theoretical SAC function Q(t) 
we adopted the analytical expression for a Gaussian line pro- 
file and an homogeneous emitting envelope given by Friedjung 
& Muratorio STW7\ . 

All the quantities are in cgs units. The term w in eq. 0] 
describes the overall deviation of the population of the upper 
(odd) levels of the Fe n transitions with respect to the popula- 
tion of the lower (even) ones. If we assume that the odd and 
even level population laws link up at an intermediate x , the 
term turns out to be equal to w - (% - Xc) (<7 _ P), where for 
the case of Fe n one can take p = 4.72 e V. Since we expect 
the odd terms to be largely underpopulated with respect to the 
even ones (see e.g. Fig. 13 in Verner et al. (I2002> '). this value 
of w should be taken as an upper limit, which corresponds to a 
lower limit to S'. 



Table 3: SAC fit parameters for the permitted Fe n narrow 
component lines, p (eV -1 ), X c , q (eV -1 ) and Y e are fitting 
parameters,^- 2 is a measure of the fitting precision. 



year 


1995 


1999 


2001 


P 


1.32 ±0.33 


1.28 ± 0.23 


0.94 ± 0.31 


x c 


-7.27 ± 0.09 


-7.52 ± 0.08 


-7.37 ±0.11 


q 


1.79 ±0.19 


1.60 ± 0.16 


1.28 ±0.17 


Y r 


-17.00 ±0.54 


-17.40 ±0.47 


-18.37 ±0.47 


N lines 


50 


50 


28 


red.^ 2 


0.65 


0.45 


0.34 



Fig-Elshows the result of the least squares fit of the dered- 
dened fluxes (with error bars) of 49 narrow components of the 
Fe ii lines measured in the 1999 spectrum of HD 45677. The 
high excitation (x u > 7.5 eV) Fe n lines have not been included 
in the fit, and will be discussed separately in Section l3~2l 

As expected, in Fig.[6]multiplet 74 emission lines lie on the 
optically thin branch of the SAC, while the very strong lines of 
multiplet 42 lines are the most optically thick ones (central line 
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Figure 6: Global SAC of the narrow emission components in 
the 1999 spectrum of HD 45677. Results from the least squares 
overlapping of multiplets 27, 28, 32, 37, 38, 41, 42, 48, 49 and 
74. The representative points are shifted horizontally and verti- 
cally according to a least squares fitting procedure by log level 
population corresponding to excitation temperatures of 3900 K 
and 3150 K for the lower and upper terms respectively. 

opacity of about 70). Note that many multiplets lie around the 
bend of the SAC (log t 0). 

The fit gives a mean lower level population p correspond- 
ing to an excitation temperature T/ of 3900^™] K, in fairly good 
agreement with the values previously found by considering two 
groups of multiplets having common upper terms. The result- 
ing mean upper level population parameter q corresponds to a 
mean excitation temperature T„ of 3150*™ K, much smaller 
than that of the lower levels. Let us note that the modelling of 
Verner et al. J2002> for two of the Weigelt blobs of 77 Car pre- 
dicted electron temperatures going down to 5000 K. 

We have applied this analysis also to the narrow Fe 11 emis- 
sion components in the 1995 and 2001-2002 spectra of HD 
45677. The derived best fit parameters for the three epochs are 
summarized in Table. [3] One may note in the table that the 
derived Y c parameters are very different from one epoch to an- 
other. This cannot be totally ascribed to a variation of the size 
of the Fe n emitting region, to which the Y c parameter is re- 
lated, since, as noted in Sec. 12.21 the spectrum did not largely 
vary during 1995-2002. Actually, the fitting value of Y c is cor- 
related with the q parameter, so that a small error on q results 
in a large change in Y c . In view of this point, as also discussed 
above, in this paper we shall mostly concern ourselves with 
the 1999 spectrum of HD 45677, and leave the problem of the 
spectral variation to a forthcoming paper. 

The same procedure has been applied to the Fe 11 broad 
emission components components. Fig. shows the global 
SAC of the Fe 11 broad emissions in the 1999 spectrum of HD 
45677 obtained from the least squares overlapping of ten Fe 11 
multiplets. The fluxes of the broad components are less pre- 
cisely determined than the narrow component ones, and the 
scatter of their SAC curves is higher, respectively. 



Table 4: SAC fit parameters for the permitted Fe 11 broad com- 
ponent lines. See Tableland text for definitions. * parameter 
taken equal to the narrow components one as no convergence 
is obtained. 



year 


1995 


1999 


2001 


P 


1.32* 


1.28* 


0.94* 


x c 


-6.97 ± 0.14 


-7.03 ±0.11 


-7.29 ± 0.18 


q 


0.99 ± 0.29 


1.50 ±0.26 


1.53 ±0.26 


Y c 


-19.35 ±0.82 


-17.84 ±0.75 


-17.84 ±0.74 


N lines 


49 


49 


28 


red.^ 2 


0.45 


0.41 


0.33 
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Figure 7: SAC of broad components of multiplets 27, 28, 32, 
37, 38, 41, 42, 48, 49 and 74 of Fe n in 1999. The representative 
points are shifted horizontally and vertically by log level popu- 
lation corresponding to excitation temperatures of respectively 
3900 K and 3400 K 



The fitting parameters of the Fe n broad emission compo- 
nents are given in Table 0] 

The SAC analysis of the broad components in 1995 and 
2001-2002 is more uncertain, due to the lower spectral resolu- 
tion in 1995, which results in a less precise determination of the 
broad component fluxes, and to the smaller set of data in 2001- 
2002. Despite the lack of convergence for the p value in 1999, 
the parameters, summarized in Table. @] appear fairly similar to 
the ones derived from the narrow emission components except 
for the 1995 data, despite the too small q value derived. 

3.2. High excitation Fe n permitted lines 

In the spectrum of HD 45677 we have identified many emission 
lines arising from high energy levels (x„ > 7.5 eV), that have 
not been included in the previous SAC analysis. Fig. [8] shows 
the position of the narrow components of these lines with re- 
spect to the SAC derived from the narrow components of the 
lower excitation Fe n lines. Indeed, plotting these lines in the 
same graph may present a difficulty, since we cannot adopt for 
these lines the same excitation temperature of the lower excita- 
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Table 5: For each high excitation line we give : (1) the labo- 
ratory wavelength, (2) upper and lower term of the transition, 
(3) upper and lower term excitation potential, (4) logarithm of 
oscillator strength, (5) log ratio of observed flux to calculated 
flux according to the lower excitation lines parameters. 



Wavelength 


upper Term 
lower Term 


Xup 
Xlow 


log gf 




6158.019 


(3G)4p x4F 
(3F)4s c4F 


8.23 
6.22 


-3.09 


4.1 


6160.755 


(3H)4p z4H 
(lF)4s c2F 


7.58 
5.57 


-3.63 


3.5 


6233.534 


d5 4s2 c4D 
(5D)4p z4F 


7.47 
5.48 


-2.94 


2.7 


6248.898 


d5 4s2 c4D 
(5D)4p z4D 


7.49 
5.51 


-2.70 


2.9 


6317.983 


d5 4s2 c4D 
(5D)4p z4D 


7.47 
5.51 


-1.99 


2.6 


6371.125 


d5 4s2 c4D 
(5D)4p z4F 


7.49 
5.55 


-3.56 


2.8 


6383.722 


d5 4s2 c4D 
(5D)4p z4D 


7.49 
5.55 


-2.27 


2.7 


6385.451 


d5 4s2 c4D 
(5D)4p z4D 


7.49 
5.55 


-2.62 


2.9 


6442.955 


d5 4s2 c4D 
(5D)4p z4F 


7.47 
5.55 


-2.88 


2.9 


6455.837 


d5 4s2 c4D 
(5D)4p z4D 


7.47 
5.55 


-2.95 


3.0 


6482.204 


(3G)4p x4G 
(3F)4s c4F 


8.13 
6.22 


-2.27 


3.3 


6491.246 


d5 4s2 c4D 
(5D)4p z4D 


7.49 
5.58 


-2.79 


3.0 


6493.035 


d5 4s2 c4D 
(5D)4p z4D 


7.49 
5.58 


-2.58 


2.9 


6506.333 


d5 4s2 c4D 
(5D)4p z4F 


7.49 
5.59 


-3.11 


3.1 


6517.018 


d5 4s2 c4D 
(5D)4p z4D 


7.49 
5.58 


-2.76 


3.1 



tion Fe n lines, or even define an excitation temperature for 
the upper levels, as they certainly are far from equilibrium. 
However, it is clear from the figure that, whatever level exci- 
tation law is adopted, the high excitation lines are 2-3 orders 
of magnitude stronger than if they were excited by the same 
mechanism as the lower excitation permitted lines. The values 
displayed in Table [5] are a clear indication of an overpopula- 
tion of the upper levels, as also appears in Fig.^] This will be 
discussed in Section l5~T1 

3.3. Fe //forbidden emission lines 

Now we analyse the narrow and broad components of the for- 
bidden lines of ionized iron observed in the 1999 spectrum of 
HD 45677. We adopted the [Fe n] transition probabilities of 
Quinet et al. (119961 . and took for simplicity Xc — eV. 

The plot of [Fe n] individual multiplets in the 
[log gf/l, log £4-] diagram shows a systematical down- 
ward shift of multiplets arising from higher excitation terms, 
that can be used to determine their upper level population 
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Figure 8: Position of the high excitation Fe n lines (filled tri- 
angles) in the 1999 spectrum of HD 45677 with respect to the 
global SAC of Fig. [6] The multiplet 37 lines (filled squares) are 
shown for comparison . 
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Figure 9: SAC diagram of narrow components of multiplets 4F, 
6F, 7F, 17F, 18F, 19F, 20F and 21F of [Fe n] in 1999. The same 
mean level excitation temperature of 10500°K has been used 
for both lower and upper levels. The observational points are 
interpolated with horizontal line. 

gradient q. In the strongest multiplets there is also an indication 
of the most intense lines having a smaller normalised flux. 

We have overlapped the [Fe n] multiplets in the 
[X - X c , Y - Y c ] plane, in the same way as we did for permit- 
ted Fe ii lines , with the least squares method, assuming that 
the ground and lower Fe + levels have the same excitation tem- 
perature (i.e., p = q = 5040/Ty). The derived SAC diagram 
shown in Fig.|5]corresponds to a best fit excitation temperature 
of the metastable levels (0-4 eV) of T f = 10500 + 1 000 K. 
Actually, the upper [Fe n] levels are the same as the lower lev- 
els of the permitted transitions treated above. Nevertheless, the 
fact that two such different excitation temperatures have been 
derived, added to the smaller FWHM value of the [Fe n], 25 
kms~' against 45 kms -1 for Fe n permitted lines, is strongly 
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Table 6: SAC fit parameters for the [Fe n] narrow components, 
q (eV~ l ) and Yj are fitting parameters. Y is Yf corrected for q 
at a mean 3 eV excitation potential, log N(H)is the logarithm 
of the total number of Hydrogen atom in the emitting volume. 



year 


1995 


1999 


2001 


q 


0.44 ± 0.05 


0.47 ± 0.05 


0.57 ± 0.08 


Yf 


-15.51 ± 0.16 


-15.41 ±0.13 


-15.20 ±0.25 


Y 


-16.83 


-16.82 


-16.91 


Tf 


1 1 500 ± 1 000 


10500 ± 1000 


9 000 ±2 000 


logoff) 


49.9 ± .2 


50.0 ± .1 


50.2 ±.3 


N lines 


85 


80 


31 


red.^ 2 


0.68 


0.82 


0.94 



Table 7: SAC fit parameters for the [Fe n] broad components. 



yc£ir 


1995 


1999 


2001 


q 


A CH i A OA 

0.50 ± 0.20 


A *C A i A AA 

U.64 ± 0.09 


(\ AO i f\ 1 O 

0.48 ±0.18 


•v 

if 


-15.95 ± 0.6U 


if to . at/; 

-15. JO ± U.ZO 


-15.94 ± U.53 


Y 


-17.45 


-17.30 


-17.38 


T f 


10 000 ±3 000 


8 000 ± 1 500 


10 500 ± 3 000 


logN(H) 


49.5 ± .6 


50.0 ± .3 


49.5 ± .5 


N lines 


61 


63 


28 


redvt 2 


2.0 


1.23 


1.45 



As for the broad components of the [Fe n] emission lines, 
the strength of the better measured lines seems to suggest by 
itself an emitting region with properties comparable to those of 
the narrow components discussed above, this is confirmed by 
the fitting values given in Table. 



3.4. Cr ii emission lines 



suggestive of two different formation regions. The results for 
the three epochs are summarised in Table[6] 

Let us point in Table [6] that the variation in the different 
epochs of the Tf and Yf parameters are much smaller than those 
derived for the permitted lines, probably as the result of the 
smaller uncertainties. The value, at the mean upper level exci- 
tation potential of our sample lines J^, — 3.0 eV, of the quantity 
Y = -p + Yf is smaller by a factor .8 in 2001 correspond- 
ing approximately to the small flux decrease between 1999 and 
2001-2002 noted in Sec. EH 
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Figure 1 1 : Self Absorption Curve for the narrow components 
of Cr n emission lines . 



Table 8: SAC fit parameters for the narrow and broad compo- 
nent lines of CrII in 1999. 



component 


narrow 


broad 


p=q 


0.82 ±0.15 


0.43 ± 0.28 


X c 


-8.63 ± 0.60 


-7.03 ± 1.00 


Y r 


-18.42 ±0.98 


-20.87 ± 1.85 


N lines 


11 


11 


red.^ 2 


0.07 


0.19 



Figure 10: SAC diagram of broad components of multiplets 4F, 
6F, 7F, 17F, 18F, 19F, 20F and 21F of [Fe n] in 1999. The 
representative points are shifted horizontally and vertically by 
log level population corresponding to an excitation temperature 
of 8 500 K. 



The analysis of the narrow and broad Cr n multiplet 44 
components, indicates a moderate optical thickness. Here again 
we adopted^ = QeV in the equations defined in Section |3~T1 
Due to the low number of Cr n lines observed in the other mul- 
tiplets we have derived the SAC parameters taking the same 
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Figure 12: Self Absorption Curve for the broad components of 
the Cr ii emission lines. 



excitation temperatures for the lower and upper terms of both 
the narrow and broad components. 

The SACs plotted in Figs. II llandll2lcorrespond to the pa- 
rameters given in Table. [8] 

4. Line profiles 

The two-component profile of the emission lines is sugges- 
tive of line formation in two regions with different velocity 
broadening in the observer's line-of-sight. In B[e] stars, of 
which HD 45677 is one important representative, low excita- 
tion emission lines, like among others Fe n and Cr e, are gen- 
erally thought to be formed in an equatorial disk around the B 
star, while the higher temperature emission lines are thought to 
be emitted in a polar wind. Zickgraf d2003t recently studied the 
kinematical structure of B[e] star envelopes based on high res- 
olution profiles of permitted and forbidden emission lines, and 
compared them with theoretical profiles emerging from an op- 
tically thin latitude-dependent expanding stellar wind and an 
equatorial opaque dust ring. He observed in his 1986-88 CAT- 
CES spectra of HD 45677 a marginal line split of 6 km s~ ! for 
the [O i] 6300 A emission line, while [Fe n] and [N n] were sin- 
gle peaked. He interpreted his observations with line formation 
in a wide opening angle polar wind seen nearly pole-on. 

In HD 45677 the profile of the emission lines do not show 
the double peaked structure typical of line emission from an 
equatorial rotating disk seen at high inclinations. The narrow 
components of strongest permitted emission lines are split by 
a slightly longward displaced ( — h3kms _1 ) dip which we at- 
tribute instead to absorption from matter flowing towards the 
star. A central absorption with the same radial velocity is also 
present in the Na i yellow doublet, and in the Mg n 448 1 A line. 
No line splitting is seen in the forbidden lines. Another impor- 
tant feature of HD 45677 is the narrower profile of the forbid- 
den [Fe n] lines, suggesting line formation in regions with a 
smaller radial velocity range, like for instance the outer parts 



of a rotating disk in which the inner, denser and faster rotating 
parts produce the permitted emission lines. 

In HD 45677 the velocity broadening of the Si n photo- 
spheric lines (<~ 70kms~ 1 , Israelian et al. I1996i should indi- 
cate either a relatively low stellar rotational velocity, or a low 
inclination angle i of a fast rotating star. If this is the case, the 
narrow emission components could arise from a dense, fast ro- 
tating disk seen nearly face-on. If one assumes that matter in 
the disk is rotating at the Keplerian velocity, the narrower [Fe 
n] emission lines should be formed in the outer (and less dense) 
rings of the disk. In such a framework, the broad components 
of the Fe n and [Fe n] emission lines, that, as discussed above, 
approximately emit a total amount of radiative power compa- 
rable to that of the narrow components, should be formed in 
region(s) with a large radial velocity range in the observer's 
line-of-sight. This could be identified as being a rather massive 
stream of matter flowing from (or to) the stellar polar regions. 

We geometrically model the emitting region of HD 45677 
in order to synthetise simultaneously both the observed profiles 
of the [Fe n] emission lines and the Fe n ones. 

As discussed above, the model of a disk seen nearly edge- 
on proposed for this star for instance by Swings [19731 or de 
Winter et al. 119971 besides the fact that it would imply a too 
small star rotational velocity, would generate emission line 
profiles with the typical double-peaked profile totally different 
from the ones observed. 

The bimodal model currently adopted for B[e] stars sup- 
poses different physical conditions in the dense disk medium, 
where the emission lines of the lower excitation species such 
as Fe ii are supposed to be produced, and in the polar wind or 
streams, where lines from higher excitation species are emitted 
(Swings 1973 Zickgraf 1998J. We nevertheless took this con- 
figuration, namely an edge-on disk and a polar wind confined 
in a cone having the same physical conditions (temperature and 
density), into consideration. The most similar profile was ob- 
tained with a 24° opening angle disk inclined 85 on the sky 
plane, in which the wide component is emitted, and a 30 cone 
wind for wich we had to suppose a 100 km s~ l decreasing ve- 
locity wind to get the narrow component. The computed profile 
is shown in Fig.^] But the result is unsatisfying (too sharp). 

In order to get better results we had to suppose that the 
narrow components are produced in a nearly pole-on opti- 
cally thin disk. We got satisfying results taking an inclination 
to the sky plane of 12°, and as a first step a keplerian rota- 
tional velocity taken to be 340 km s~ l at the star surface. With 
such an inclination and velocity we are in agreement with the 
v sin i = 70 km s~ l determination of Israelian et al. J1996i . 
Such a region of line formation, would, when seen nearly pole 
on, resemble Kastner's static model. 

To get the broad line components we had to suppose the ex- 
istence of a flow falling on the disk (or alternatively escaping 
from the disk), whose material is incorporated into (or alterna- 
tively ejected from) the disk material at the Keplerian velocity 
in each point. 

With such a configuration we can reproduce at the same 
time as well the permitted (Fig. 1141 as the forbidden ("Fig. I15i 
Fe ii line profiles, by just playing on the distance to the star of 
the emitting zone. 
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Figure 13: Histogram of the emitting material radial velocities 
of an optically thin disk plus polar wind. The disk is inclined 
85° to the sky plane with a total opening angle of 24°. The ve- 
locity in the disk is 150 kms~' at the photospheric radius r=l, 
and the law is Keplerian. The polar wind is emitted in a 30° 
cone with a velocity of 100 kms~' at the pole and decreases 
exponentialy. In this figure and in figures ^] and the syn- 
thetised profile and the continuum are plotted as a solid line, 
each Gaussian component as a short dashed line, the total fit 
profile as a dotted line and the residual flux as a long dashed 
line. 



The permitted and forbidden FWHM line components can 
be understood if we assume that the permitted lines are formed 
in the inner parts of the disk, while the forbiden ones issue from 
the external parts of the Keplerian disk. 

Fig. [21 shows the profile corresponding to an optically thin 
Keplerian disk inclined 12° on the sky plane. The velocity at the 
photospheric radius r=l is equal to 340 kms~'. The emission 
comes from the inner part (5 < r < 9, in photospheric radius 
unit) of the optically thin disk whose total opening angle is 
24°. The components flux ratio narrow/broad is close to 1 (the 
narrow component is slightly more intense) as was found for 
HD 45677, while the components FWHM are similar (44 and 
194 km s ) to the Fe n permitted lines ones. 

Fig. 1151 shows the profile corresponding to the same disk, 
but with emission coming from outer regions 
(15 < r < 20, in photospheric radius unit). The component flux 
ratio narrow/broad is slightly higher than in Fig.^] as it is the 
case for [Fe n] with regard to Fe n lines in HD 45677, while 
the components FWHM are similar (28 and 108 km s _1 ) to the 
[Fe n] lines ones. 

These two profiles model fairly well the mean component 
parameters we measured in HD 45677 as well for the permitted 
as for the forbidden Fe n lines. 




-100 100 



Figure 14: Theoretical profile corresponding to the histogram 
of the emitting material radial velocities of an optically thin 
Keplerian disk inclined 12° on the sky plane. An accreting 
wind generates the broad component. The total fit profile co- 
incides with the broad component for large velocities, and with 
the sum of the components for small velocities. The component 
FWHM are the ones derived in Sec. 12.21 The Keplerian veloc- 
ity is equal to 340 kms 1 at the photospheric radius r=l. The 
emission comes from 5 < r < 9, in photospheric radius unit. 
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Figure 15: Theoretical profile corresponding to the histogram 
of the emitting material radial velocities of an optically thin 
Keplerian disk inclined 12° on the sky plane. An accreting wind 
generates the broad component. The component FWHM are 
the ones derived in Sec. 12.21 The Keplerian velocity is equal to 
340 km s _1 at the photospheric radius r= 1 . The emission comes 
from 15 < r < 20, in photospheric radius unit. 
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Figure 16: Relative population of the Fe + levels in the [Fell] 
(squares), and in the low and high excitation Fe II emission re- 
gions derived from the narrow emission components measured 
in the 1999 spectrum of HD 45677. The fluxes of the lower ex- 
citation emission lines (open circles) have been corrected for 
self-absorption, as discussed in the text. The high excitation 
emission lines are drawn as triangles. The straight lines rep- 
resent the mean excitation temperatures of the upper levels of 
the [Fe n] and of the lower and upper levels of the Fe n emis- 
sion lines, respectively. The crossing at 4.72 eV marks the con- 
nection between the lower (3900 K) and upper (3150 K) Fe n 
level population laws. The zero points Y(0 eV) are the esti- 
mated ground level relative populations. 

5. Discussion 

5.1. Level Population 

The results of our analysis concerning level population are 
summarised in Fig. 1 161 where we plot for the three groups of 
lines - forbidden, low excitation and high excitation permit- 
ted Fe ii narrow components emission lines - the logarithm of 
the normalised line flux as a function of the excitation poten- 
tial Xu of the upper level of the transitions for our 1999 data. 
The fluxes of the intermediate excitation permitted lines are 
corrected for the self-absorption as discussed above. The di- 
agram is used to investigate the formation region of the three 
line groups. One normally expects the upper levels of the per- 
mitted lines to be de-populated with respect to the metastable 
levels. If no de-population occurs, the population laws for the 
upper levels and the lower levels should join. In fig. ^] the 
mean population laws for the upper (3 150 K) and the lower 
(3 900 K) cross at the mean excitation potential value of 4.72 
eV, to take into account the overall deviation of the population 
of the upper (odd) levels of the Fe n transitions with respect to 
the population of the lower (even) ones (see Section l3~T1 l. The 
extrapolated ordinate value at eV of this law, given by equa- 
tion 0] for the Fe n lines (-14.51 ± 1.2) turns out to be close 
to the value for the [Fe ii] lines (-15.41 + 0.13). This will be 
discussed in Section l5~2l 

We get in this diagram a limit to the position of the 
metastable levels in the permitted line formation region. As for- 



bidden and low excitation Fe n transitions have the metastable 
levels in common, we may note that the [Fe n] transition upper 
levels in the Fe n permitted lines region appear to lie in this di- 
agram aproximately 2 to 3 dex under the location of the [Fe n] 
observed line ones. Consequently [Fe n] lines in the permitted 
line formation region cannot be observed. 

Finally, the high excitation Fe + levels are clearly overpopu- 
lated with respect to the upper levels of the lower excitation 
Fe ii emission lines. We attribute this overpopulation to the 
same fluorescence mechanism acting in other peculiar stars, 
such as KQ Pup, 77 Car, and RR Tel (Muratorio et al. WWh 
Viotti et al. 119991 Hartman & Johansson 2000) whose optical 
and UV spectrum shows prominent Fe 11 lines arising from high 
excitation levels. Like in these cases, the far-UV radiation is 
most probably an important process of population of the Fe + 
upper levels in the envelope of HD 45677, while probably di- 
electronic recombinations are unimportant. 

5.2. Column density 

From the position of the bend in the SAC plots of Fe 11 lines 
(Fig. |6] and Fig. □ for 1999 data) we are able to determine 
the value of X c and hence through eq. [3] the values of the 
Fe + column density of the Fe 11 permitted lines emitting regions 
for each studied epoch. 

To derive the Fe II column density ^ , defined by eq. (3), 
in the narrow component Fe 11 permitted line region, we have 
to estimate the velocity broadening of the opacity profile v c 
in eq0 The narrow components of the emission lines of HD 
45677 should be formed in a region with a small velocity gra- 
dient in the line of sight, as discussed in Section |4] it is sup- 
posed to be an equatorial disk seen nearly pole on. If we are 
observing a disk in rotation, without being in its plane, the ap- 
propriate v c to use, to which the line broadening along a line 
of sight is linked, is twice the thermal broadening velocity that 
has been derived from the [Fe 11] emission lines, as discussed 
in Section 13.31 The local velocity dispersion for a tempera- 
ture T th = 1 1 500 K is given by v th = 0.2 1 5 1 5 VWPA where 
PA = 55.85 for Fe + (Baratta et al. 119981 . For the narrow com- 
ponents v c = 6.2 km s~ l . As for the broad components, the ap- 
propriate v c to use is given by our model, as the mean velocity 
along a line of sight, v c is equal to 7 km s . 



Table 9: logarithm of the minimum Fe + column density -j- (for 
an excitation potential^ = 2.75 eV) of both components emis- 
sion line region at the different epochs studied. The first line 
display (N) the narrow component region values, while the sec- 
ond line displays (B) the broad component region ones. 



year 


1995 


1999 


2001 


log £ (N) 


14.6 + 0.1 


14.9 ±0.1 


14.73 ±0.1 


log ^ (B) 


14.4 ±0.1 


14.5 ±0.1 


14.5 ±0.1 



The different epochs and region log y- values at Xc — 
2.75 eV are given in Table|9] If we extrapolate this column den- 
sity to eV using 1.28 as the exponent of the power law of the 
lower terms population we obtain log N(Fe + ) = 46.5 in 1999. 
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It is not certain that assumption (d), in sec. 3, about the same 
Boltzmann population distribution for all metastable levels is 
correct for the more excited of these levels (see Verner et al. 
I2002> . Therefore, if we are very conservative and suppose that 
the mean excitation temperature down to ev is 10 000 K, typi- 
cal of HII regions, as well as being of the order of what we find 
here for the less excited forbidden [Fe n] lines with lower levels 
at or near ev, the column density values might be smaller by 
a factor of about 2.1 dex. 

The different narrow and broad component column densi- 
ties in Tableware very similar. The eventual column density 
variations with time are of the order of the precision. 

As for forbidden lines, Fig. [9] shows a marginally signif- 
icant bending of the data points towards large optical thick- 
nesses, as if the strongest [Fe n] lines were slightly self- 
absorbed. The curve bend, if real, would correspond in this 
figure to a self absorption factor of 3 for strongest of the ob- 
served forbidden lines, much less than what is observed for the 
strongest permitted lines, as the optical depth range is lower for 
forbidden lines. It would correspond to an Fe + column density 
about 100 times larger than that above derived for the permitted 
lines, which seems to us unrealistic. 

The vertical overlapping of the multiplets provides the zero 
point of the ordinates: Yf = -15.41 + 0.13 (see Table[6j. This 
quantity is related to the total number A^(Fe + ) of Fe + ions in the 
[Fe n] emitting volume through the relationship (see Baratta et 
al. 1998): 

log jV(Fe + )/d 2 = Y f + log U+ 16.977 (5) 

, where U is the partition function of Fe + . 

Assuming for HD 45677 a distance of 420 pc (Zorec [T998l . 
we obtain log A^(Fe + ) = 45.6 ± 0.13 which gives in 1999, for a 
cosmic abundance, a total hydrogen mass of 1.7 10 26 g, with a 
precision factor estimated as 0. 13. 

Notice that the logA^(Fe + ) = 46.5 ±1.3 value of the Fe n 
permitted emitting region, obtained with a value of 1.28 for the 
lower level constant p of eq. (1), is close to the forbidden Fe n 
region one. 

5.3. Size of the emitting region 

From the fit results summarised in tables[3]and|4]we can calcu- 
late, using equations[3]and|3] the apparent size of both emitting 
line regions. 



Table 10: Emission line region minimum sizes at the different 
epochs studied (in cm). The first line display the narrow com- 
ponent region values (N), while the second line displays the 
broad component region ones (B). 



year 


1995 


1999 


2001 


logR'(N) 


13.1 ± 0.5 


12.9 ±0.3 


12.5 ±0.4 


logR' (B) 


12.8 ±0.8 


13.0 ±0.4 


12.6 ±0.5 



For 1999 we obtain logR' = 12.90 + 0.34 for the narrow 
line region. The values of logR' in Table. [TUI are typical radii 



of the line emitting region, equal to J Sc ™' , The radius of 
a B2 V star is according to Drilling and Landolt (2000 in 
'Astrophysical Quantities", page 389)) near 4 10 11 cm. The size 
of the region where the narrow Fe II emisison lines originate, 
probably a disk, turns out to be 20 (from 9 to 44) times that of 
the stellar radius. The model we presented in Section[4]is con- 
sistent with these results as the permitted line profile similar to 
the one observed is emitted from a region whose surface is that 
of a disk of 7 stellar radii. The assumption has been made in the 
R' determination that the population laws for the upper levels 
and the lower levels join at 4.72 eV. Their values obtained are 
therefore minima, corresponding to maximum w values. 

The logR' values obtained from 1995 and 2001 data in 
Table [TUI remain close to the limits evaluated from the 1999 
data in Table ITU1 

As for the broad component emitting region apparent sur- 
face, the values of log R' are very close, at least in 1999 to the 
narrow component one. This result is consistent with the model 
of a disk seen nearly pole-on, developped in Section [4] as the 
regions in which each component is emitted, being seen under 
the same opening angle, are supposed to have similar apparent 
surfaces. 

As for the narrow component Cr n emitting region, we used 
v'o = 6kms~ l (twice the 10 000 K thermal broadening velocity) 
and the values of Table. [3] hence a CrII column density — = 
15.98 + 0.6. The logR' value derived in 1999 (12.2 + 0.8) is 
only slightly smaller than that of the Fe n emitting regions. 

6. Conclusions 

The SAC method has given consistent results. We have ob- 
tained the properties of what is probably a disk of HD 45677, 
including limits on a characteristic column density and a char- 
acteristic radius of the Fe n emitting region. These estimates re- 
quire reasonable assumptions, one of which is a true Boltzmann 
distribution of the populations of the lower even metastable lev- 
els and another that the Boltzmann-like law of the populations 
of the odd upper levels of strong optical lines are not more than 
what would be expected from a smooth fit between these odd 
and even population laws. We have demonstrated that the well 
known split emission line profile is in fact due to absorption 
components like in the ultraviolet. We have demonstrated that 
emission lines have narrow and broad components. We have 
shown that the double emission profile can be explained with 
a disk plus a flow (perhaps a wind) moving in or out of the 
disk. The Fe n forbidden lines show a similar two component 
structure, but the smaller width of each component is consistent 
with their production in a more external part of the disk, where 
the Keplerian velocity is lower. The examination of the fluxes 
of the narrow [Fe ii] components has shown that the metastable 
Fe + levels follow a Boltzmann-type law with a mean excita- 
tion temperature around 1 1 000 K. Such a value may actually 
represent the electron temperature of the [Fe n] emitting re- 
gion, if the metastable levels are in equilibrium with the ground 
term through electron collisions. This is probably the case of 
HD 45677, since we expect the electron density to be large 
enough (N c > 10 6 cm -3 ) for bringing the metastable levels 
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to LTE (Nussbaumer & Storey 1988). The low electron tem- 
perature with respect to the stellar surface temperature may be 
due to cooling of the [Fe n] region by the large line emission 
of this ion. Indeed, low electron temperatures were predicted 
by Verner et al. J1999I and Verner et al. J2002> for the 77 Car 
Weigelt blobs that are strong Fe + emitters. Fairly low [Fe 11] ex- 
citation temperatures were derived for 77 Car (12000 K, Viotti 
et al. JT9991 and for the symbiotic star RR Tel (6600 K Kotnik- 
Karuza et al. l2002l . 

The narrow components of the Fe 11 permitted emission 
lines in HD 45677 are slightly, but significantly broader than 
the corresponding [Fe 11] components. From the SAC analysis 
of the Fe 11 emission lines we have derived an Fe + lower level 
excitation temperature much smaller than that derived from the 
[Fe 11] lines, namely around 3 900 K. We take these results as 
evidence of different (but possibly partially overlapping) emit- 
ting regions. The derived excitation temperatures of the lower 
and upper levels of the permitted Fe 11 lines are much smaller 
than that expected if Fe + were in thermodynamical equilib- 
rium with the local plasma. Verner et al. (|2002f calculated the 
formation of singly ionized iron in the Weigelt blobs of 77 Car. 
According to them Fe + emission is the most important coolant 
in the blobs. Fe 11 emission originates at electron temperatures 
5 000 K < T e < 7 500 K. This can also be seen by comparing 
their figs, 8 and 9. We argue that for HD 45677 the upper lev- 
els of the permitted low excitation transitions are mostly pop- 
ulated by the diluted stellar UV radiation followed by sponta- 
neous decay, and that in the Fe 11 region the radiation density is 
high enough to make the radiation induced transitions more im- 
portant than electron collisions. Similar considerations can be 
made for Fe 11 and [Fe 11] broad components, as their emitting 
region parameters appear fairly close to the narrow components 
ones (see Tableland Table. 0]for Fe 11 lines and Table. [6]and 
Table. 0for [Fen] lines). 

The Fe 11 metastable level, as well as the lower level 
(Table. [5} excitation temperatures are suspected to increase 
with time. The high excitation anomalous multiplets may in- 
dicate the presence of inhomogeneities. 

The shape of the iron broad component is insufficiently 
known now to decide about the structure of the non-disk com- 
ponent. Nevertheless the similarity of the physical parame- 
ters derived for narrow and broad components (Texc, log -^) 
suggest that both emitting regions, if not the same (different 
FWHM) are physically linked. 

Our identification of the broad components, and the sus- 
pected variability of both emitting regions physical parameters 
should stimulate in future new higher quality observations (es- 
pecially in the blue) in order to study more precisely the shape 
of the emission lines. 
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